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Abstract 

We investigate 'Scaling ansatz' in the neutrino sector within the framework of type I seesaw 
mechanism with diagonal charged lepton and right handed Majorana neutrino mass matrices (Mr). 
We also assume four zero texture of Dirac neutrino mass matrices (mo) which severely constrain 
the phenomenological outcomes of such scheme. Scaling ansatz and the present neutrino data allow 
only Six such matrices out of 126 four zero Yukawa matrices. In this scheme, in order to generate 
large #13 we break scaling ansatz in through a perturbation parameter and we also show our 
breaking scheme is radiatively stable. We further investigate CP violation and baryogenesis via 
leptogenesis in those surviving textures. 

PACS numbers: 14.60.Pq, ll.30.Hv, 98.80.Cq 



I. INTRODUCTION 



Neutrino physics is now playing a pivotal role to probe physics beyond the Standard Model. 
Confirmation of tiny neutrino masses as well as nonzero mixing angles have thrown light on 
the structure of the leptonic sector. In the quest towards understanding of a viable texture of 
neutrino mass matrix popular paradigm is to advocate flavor symmetries, directly associated 
with some gauge group. On the other hand, there are some other ansatzs which also give 
rise to interesting phenomenological consequences, although their origin from a symmetry 
discrete or continuous are yet to be established at the present moment. 

In the present work we bring together two ideas to explore the neutrino phenomenology, 
particularly, to generate large #i 3 [Tj- [T8] as reported by recent experiments [19]- [21] as well 
as CP violation and baryogenesis via leptogenesis. In this scheme we consider 

a) Scaling ansatz |2"5]-[3"3"]. 

b) Four zero texture [31]- [38] of Dirac neutrino matrix (mc), 
within the framework of type I seesaw mechanism denoted as 

m v = -m D M- R l m T D (1) 

where Mr is a 3 x 3 right chiral neutrino mass matrix and we consider the basis in which 
charged lepton and Mr are flavor diagonal. 

Scaling ansatz [25]- [33] posses a distinctive feature that the texture is invariant under renor- 
malisation group evolution unlike other symmetries such as fi — r symmetry. Basically, the 
ansatz correlates the elements of neutrino mass matrix through a scale factor and it can be 
implemented in different ways. Although the theoretical origin of such ansatz is not yet well 
known, however, this ansatz can be approximated as S2L symmetry (i.e fi — r symmetry in 
the left handed neutrinos) with the value of the scale factor unity. Furthermore, it leads to 
inverted hierarchy of neutrino mass with 777,3 — 0, and # 13 = 0. Thus, it is obvious to break 
such ansatz in order to generate nonzero #13. 

The other assumption that occurrence of four zeroes in mo gives rise to a more constrained 
feature that the phases contributing to the high scale CP violation required for leptogenesis 
(basically the phases of momj, matrix) are determined in terms of the low energy CP 
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violating phases ( i.e phases of m u ). 

We divide all the 9^ 4 = 126 four zero textures in the following Classes : 

i) det(m£>)= and no generation decouples: 27 textures 

ii) det(m£>) ^ and no generation decouples: 72 textures. 

iii) det(m/))= and one generation decouples: 18 textures 

iv) det(m£>) ^ and one generation decouples: 9 textures 

The textures belong to Class (ii) are already studied extensively [34j-[38j. Class (iii) and 
(iv) are incompatible with the neutrino experimental result. The remaining Class, Class 
(i), which is yet to be explored, posses one zero eigenvalue which is still allowed by the 
present experiments. The interesting point is to note that if we insert scaling ansatz to all 
four zero textures and consider those textures in which four zero remain four zero and no 
generation decouples, we see that the survived textures are only from Class (i). Motivated 
with this unique selection property of scaling ansatz, in the present work we investigate 
textures belong to Class (i). In addition to one eigenvalue zero, scaling ansatz also dictates 
one mixing angle to be zero. We further generate nonzero # 13 through the breaking of 
scaling ansatz due to a small perturbation parameter in mp. We investigate all possible 
cases and finally we demonstrate that the broken scaling ansatz textures remain invariant 
under renormalization group (RG) evolution. 

Our plan of this paper is as follows : In Section II we discuss different types of scaling ansatz 
and allowed four zero textures. Section III contains parametrization and diagonalisation of 
neutrino mass matrix. Breaking of scaling ansatz and generation of nonzero #13 are discussed 
in Section IV. Numerical results are given in Section V and Section VI contains the possible 
baryogenesis via leptogenesis scenario arises in those textures and summary of the present 
work is given in Section VII. Discussion on RG effect is given in Apendix A and explicit 
expressions arise in Section IV are included in Appendix B. 
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II. FOUR ZERO YUKAWA TEXTURES AND SCALING ANSATZ 



A. Scaling ansatz 

Several authors [25]- [33] have been studied scaling ansatz through its implementation 
along the columns of effective m u matrix. In the present work, we consider this ansatz at a 
more fundamental level of |25j and we find that implementation of this ansatz along the 
rows of mu with a diagonal Mr effectively gives rise to the same structure of m v [26] after 
invoking type-I seesaw mechanism. According to this ansatz elements of a row (of 3x3 m^) 
are connected with the elements of another row through a definite scale factor. In case of 
3x3 mo there are three types of this ansatz which are given as follows: 
i) Second and third row are related through a complex scale factor k as 



m D 



fit 



krriD 1 



(2) 



where i is column index, i = 1, 2, 3. Invoking type I seesaw mechanism 

( m ») m = -( m D)^M^m T Dja 



T 

Dja 



= -k(m D ) rj M R l j m 
with a = e, fi, r we obtain the following scaling relations in m h 
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(3) 



(4) 



We discard the other two cases where the scale factor relates ii) First and third row and iii) 
First and second row because in those cases either # 12 or 623 is zero at the leading order. 

B. Four zero Yukawa textures 



We start with a general scaling ansatz invariant mo matrix on which we will assume 
four zeroes and explore all the possibilities. Explicit structure of mc according to eqn.(j2]) 
is given by 



( di <2 2 03 ^ 
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TABLE I: Four zero Yukawa textures compatible with Scaling ansatz 



We categorise all possible four zero textures compatible with Scaling ansatz in three different 
cases as shown in Table I. The following points to be noted : 

1. We find that out of 126 four zero textures, imposition of scaling ansatz reduces dras- 
tically the number to only 12. 

2. We ignore Category B because it is not possible to break scaling ansatz keeping the 
pattern of mo matrices unaltered. Let us assume the breaking is incorporated as 
k k(l + e), the structure of all mo remain same and still invariant under scaling 
ansatz. Thus, to break scaling ansatz in Category B, we have to have reduce the 
number of zeroes which is beyond our proposition. 
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3. We also discard all the textures in Category C since one generation is completely 
decoupled from the other two which give rise to two mixing angles zero. 

Hence, the number of surviving texture is only six and all of them are from Class (i) described 
previously in the Section I. 

For Category A as the second and third row of the matrices are connected through a scale 
factor, from now on we express them as follows 



( a ^ 



kb kc 
6c 



m 11A - 



( a ^ 

kb kc 
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( a ^ 



J II A 



kb kc 
6 0c 



(6) 



0a 
kb kc 
6 0c 



\ 



m VA 



a 
kb kc 
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m VIA 



a 
kb kc 
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where a, 6, c and k are all complex parameters. 



III. PARAMETRIZATION AND DIAGONALISATION 



Parametrization 



(7) 



We parametrize the m v matrix arises after seesaw for all mo matrices in Category A in a 
generic way as 

( 1 kpe ie pe id ^ 

m v = m kpe ie k 2 (q 2 e 2 ^ + p 2 e 2W ) k(q 2 e 2 ^ + p 2 e 2W ) (8) 
v pe w k(q 2 e 2 ^ + p 2 e 2W ) q 2 e 2i ? + p 2 e 2id } 
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with the definitions of the parameters for six consecutive cases as 
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Mi c 
M 3 a 



M 3 b 
Mi a 



Mi c 
W 2 a 



M 2 b 
Mi a' 

Considering complex k as ke 10 and to as moe 10 ™, we rotate the matrix m v by e~ %[ 
diag(l, e~ % ( e+6k \ e~ l6 ) from both sides and get the m v free from redundant phases as 

^ 1 kp P 

kp k 2 re ta kre ia 
p kre ia re 1 



(9) 

% l 2 x 



TO,. 



to 



V 



(10) 



/ 



where 



q 2 e m-e) +p2 = reia 



(11) 



Here to , k, p, r all are real positive parameters. We construct the matrix h(— m v mj) to 
calculate the mixing angles and mass eigenvalues. Expression of h is obtained as 

( A k\B\e-^ |S|e - ^ 
\B\e^ k 2 C kC 
\B\e^ kC C 



h = m v rrv v = TOq 



V 1 



(12) 



/ 



where 



A = 1 + k 2 p 2 + p 2 

B = \B\e i(p =p + k 2 pre ia + pre { 

C = p 2 + k 2 r 2 + r 2 

r sinafl + k 2 ) 

tan = 



1 + r cos a(l + k 2 ) 



2\ ■ 



(13) 



7 



Again factoring out the phase in h as h — > diag(e lr ^, 1, 1) h diag(e 1, 1), finally, we obtain 



h = m^A 



( 1 k\B'\ \B'\\ 
k\B'\ k 2 C kC 
\B'\ kC C 



where \B'\ = ^ and C — ( . 



il c 



(14) 



B. Diagonalization 



Diagonalizing the matrix h given in eqn.(14) as U^hU = diag(m 2 , m|, m|) we get 



m 2 A(- 



Pi - y/Pf - 4gi 
2 

P 1 + V / P 2 -AQ 1 



(15) 



where 

Pi = 1 + C"(A; 2 + 1), Qj = (A; 2 + 1)(C" - \B'\ 
and the mixing matrix is 



U 



( 



\ 



Cl2 Sijj 

~-5i2C23 C12C23 S23 

S12S23 -C12S23 c 23 



where c« = cos6L- and = sin^,-. The three mixing angles are 



tan 



2:i 



1 

1 



tan# 
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2|P / | V / TTP 

C'fl + A; 2 ) - 1 



#13 = 



(16) 



(17) 



(18) 



and the mass squared differences are 



Am 2 21 = m 2 AJP 2 - AQ 



Am 2 32 



-m 2 A( 



Pi + y/P? - 4Qj 



(19) 
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FIG. 1: Allowed parameter space for ansatz conserving four zero cases given in Category A 

In FigjTJ we plot the parameter space varying another model parameter a within the range 
—7i < a < 7i satisfying the following 3cr experimental ranges of neutrino data [39THT] 



35.5° < 6 23 < 53.5° 
31.7° < 0i2 < 37.7° 

6.90 x 10- 5 eV 2 < (Arr4) < 8.20 x 10~ 5 e\/ 2 
-2.73 x 10- 3 eV 2 < (Am 2 2 ) < -1.99 x lO^eV 2 . 



(20) 



We have also used cosmological bound on the sum of the neutrino masses as Emj < 0.5eV[321- 
Hi], and the lower bound obtained from neutrinoless double beta decay (/3/3q v ) as < 
0.35eV| 



IV. BREAKING OF SCALING ANSATZ AND GENERATION OF NONZERO 9 

We want to break the scaling ansatz in such a way that 

• becomes nonzero. 

• Four zero structure is also retained. 



The second assumption rules out all Category B textures as we have mentioned ear- 
lier. Breaking of scaling ansatz can only be incorporated in the remaining six four zero 
textures in Category A and after breaking the scaling ansatz by a dimensionless real 
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parameter e their structure come out as follows 



/ 



m 



I A 
D 



a 
kb(l + e) kc 
b c 



\ 



/ 



II A 



,m D 



a 

kb kc(l + e) 
Ob c 



( 



,< IA 



J 



( 



m 



IV A 
D 



a 
kb kc(l + e) 
6 c 



VA 



a 
kb(l + e) kc 
b c 



KL4 



a 
fc6(l + e) A;c 
b c 

a 
£;& fec(l + e) 
6 c 



(21) 



Theses structures of are free from RG effects which we have discussed in Appendix-A. 
Moreover, the breaking considered here are the most general which can be understood as 
follows: Consider the matrix m 1 ^ in which the breaking scheme is incorporated as 



( m o) M2 = k(l + e){m D ) 



t2 



(22) 



while 



( m o) M 3 = k ( m D) r3 - (23) 

Now, redefining the parameters k(l + e) — > k and — e — > e it is equivalent to break the ansatz 
in (wid)^ and (j7i£)) T3 elements. Proof of this equivalence is similar for other remaining five 
mo matrices. 

The effective neutrino mass matrix m v is same for all of them and is given by 



/ 



1 kp + kpe p 

kp + kpe k 2 re ta + 2k 2 p 2 e kre ia + kp 2 e 



\ 



V 



re 



(24) 



/ 



p kre ia + kp 2 e 

with the same definitions of the parameters (k, p, r, a) that we have already used in eqns.Q 



and (11) 



We now rewrite this m v by breaking it in two parts, one e dependent and the other inde- 
pendent of e, i.e 



' 1 kp p ^ 



kp k 2 re ia kre ta 
p kre ia re ia 



i kp » 



+ em 



kp 2k 2 p 2 kp 2 
kp 2 



(25) 
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where we have denoted the first matrix in the right hand side of the above equation by 
and the second one by m! v . Computing h t using the above m u , we get 



h t = m v m\, = m®m® + e{mP u m'J + m' v m°J) = h° + eh v 



(26) 



neglecting 0(e 2 ) terms. It is to be noted that h° is same as h, that we have obtained in 
eq.(12). After rotating out the phase appearing in h° we are left with 

/ a 1.1 n I i di\ 



h' 



A k\B\ \B i 
k\B\ k 2 C kC 
\B\ kC C 



eh" 



(27) 



where h" = diag(e ic ^, 1, 1) h p diagie^, 1, 1) and h' t = diag(e i(t> , 1, 1) h t diagie'^, 1, 1). To 
diagonalise h' t we first rotate this matrix with unperturbed diagonalising matrix U in eq. 



(17) with angles in eq. (18). The first part of h' t becomes diagonal, however, the h" part is 
not. Performing the operation U^h' t U we get 



/La o o 



h t " = rfhLU 



m 2 2 




+ e 



y* w q 
y z * q* 0J 



(28) 



where different elements of the the 2nd matrix are obtained from the explicit multiplication 
U^h^U. To diagonalise the second matrix of h" we further require the matrix 



U e 



A 



(29) 



1 ea eb 
—ea* 1 ec 
-eb* -ec* 1 J 

Explicit expressions of parameters x, y, z, q and w are given in Appendix B. We demand 



V 



that upto 0(e)the above matrix diagonalises h" of eq.(28), i.e after the operation Ujht"U £ 
the off-diagonal elements of the resulting matrix are zero and solving those equations we 
find out the unknown variables a, b, c. They come out as 

y 



{m\ — mf) 

z 

mi 



c = — - 



mi 



(30) 
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As a result of this rotation by the matrix U e we get 



/ 2 

' mi + ex 



u\h'lu e 



\ 





ml + ew 




(31) 



In a concise way, we actually have done the following 
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(32) 



where m' 2 , m' 3 are the new mass eigenvalues. m' 3 is still zero even after breaking of scaling 
ansatz because one column remain zero for all allowed mo- Hence, the total diagonalisation 
matrix in our scheme is V — UU e . Explicitly V is given by 

-C23S12 + C 12 C2?,{ m l y * mi ) -C 2 3Sl2( m 2? m 2 ) + C 12 C 2 3 C 23 S 12 (^) - Ci 2 C 23 (^ 



S23«12 - S23C121 



S 2 3S 12 ( 



+C23C 



t 2 — '"1 
+<*s(£ 



S23C12 -S23«12(^) + «23C12( 
+C23 



(33) 



To find out the three mixing angles we have to compare V with PMNS matrix. The £/pmns 
is given by 
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" s 12 c 23 



(~J rJ rJ 

c 12 6 23 6 13 e 
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~12 C 23 S 13 C 
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(34) 



(with = cos 6[j, = sinO'ij, 5 is the Dirac phase and cum, Pm are the Majorana phases.) 
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After neglecting the higher order terms in e the modified mixing angles are given by 
tan 23 = W\ ~ *23 + e(l + 4)(^§i?e^) - ^|i2e(g)) 



V33I " mf m 

tan^ 12 = §4-^ + e ( 1 + ^)^ e(2/I 



|y u | iZ v 12J {ml-m\) 



sin ^ 3 = |V 13 | « eJ(^Re(z) + ^|i2e(g)) 2 + (^|/m(^) + ^/m(g)) 2 
w 777.7 mo mf w> 



(35) 



and the CP violating phase S is given by 



^lm(z) + ^lm(q) 

tano = , , : ^ , , 36 

^Re(z) + %Re(q) 

where Cy = cosfty, Sy = sin 6^-, = tan^- are the mixing angles and m 1; m 2 are the 



masses for the ansatz conserving case of eqn. (18) and eqn. (15) respectively. From eqn. (31 ) 
we have the mass squared differences: 

(Am 2 i)' = Airij! + e(w — x) 

(AmlJ = Am\ 2 - ew (37) 
where Am^ and Am^ are mass squared differences for unperturbed scaling ansatz as in eq. 



(19). The measure of CP violation is understood through Jqp which is defined as 



j (h t ) 12 (h t ) 2 3(ht) 3 i 

CP (Ami)' (Aml 2 )' (Ami)' [ ] 

which is known function of k, p, r, a and e. 

V. DISCUSSION OF NUMERICAL RESULTS 

We explore the parameter space of the above case using the same 3a values of neutrino 



experimental data given in eqn. (20). The Lagrangian parameters p and r are ranging from 
zero to some positive values since we have separated out the phase part from them. The 
scale factor k should not have zero value because in this case the second row of (mo) is zero 
which in turn decouple the second generation. The constrained parameter space we obtain 
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FIG. 2: Allowed parameter space for e = 0.1 

as 

< p < 4 
0.65 < k < 1.4 
< r < 0.7 

-180° < a < 180°. (39) 

The values outside this range is not admissible within the above mentioned experimental 
ranges. It is to be noted that allowed parameter space in k-r plane for the ansatz breaking 
case is much larger than that in the ansatz conserving case. First of all, we found that 
throughout the allowed parameter space Em; and m v pp are always far below the experimen- 
tal bounds which could be hardly tested in the near future experiments. Next, it is amply 
clear from the expression of #13 that it is directly proportional to the value of the ansatz 
breaking parameter e. The parameter e is varied upto a reasonable choice e = 0.1 for which 
a large 6*13 is generated, however, for a smaller value of e such as e ~ 0.07, #13 ~ 10° is also 
admitted because present experimental bound on 9 13 is 3.75° < 9 13 < 13.60° for 3a bound 
from RENO [21] and 4.90° < 9 13 < 11.51° for 3a bound from Daya-Bay [23J . We have 
shown all plots for a representative value of e = 0.1 . The allowed Lagrangian parameter 
space is plotted in Fig|2j From Figj3]it is clear that #13 is almost insensitive to #23, however, 
significantly related to the values of 9 12 which is depicted in Fig.3. The CP violation param- 
arises due to nonzero #13 is plotted with #13 in Fig|4j Sign of S is not constrained 
from oscillation experiments it needs further calculation of baryon asymmetry. Plot (FigjlJ) 
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FIG. 3: Allowed values of the mixing angles for e = 0.1 
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FIG. 4: Allowed | J C p\ vs 6 13 (left), \6\ vs 0i 3 (right) for e = 0.1 

of ^13 vs 5 shows that S is maximum for smaller values of 9 13 and for larger values of 9\s, 
5 is relatively small. If we restrict #13 in 3a experimental range we have the bound on S, 
< 5 < 35° and J CP , < J C p < 0.02. 



VI. LEPTOGENESIS WITH BROKEN SCALING ANSATZ 
A. General discussion on Leptogenesis and Baryogenesis 

Let us briefly discuss about right handed Majorana neutrino decay generated leptoge- 
nesis. There is a Dirac type Yukawa interaction of right handed neutrino (iVj) with SM 
lepton doublet and Higgs doublet. At the energy scale where SU{2)l x U(1)y symmetry is 
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preserved, physical right handed neutrino Ni with definite mass can decay both to charged 
lepton with charged scalar and light neutrino with neutral scalar. Due to the Majorana 
character of Ni, conjugate process is also possible. If out of equilibrium decay of N in 
conjugate process occur at different rate from actual process, net lepton number will be 
generated. The CP asymmetry of decay is characterized by a parameter £j which is defined 
as 

We are working in a basis where right handed neutrinos have definite mass as , M B = 
diag(Mi, M 2 , M 3 ). Now, the decay asymmetry for N decay occurs at one loop level. 
Interference of tree level, one loop vertex and self energy diagrams generate the following 
for hierarchical right handed neutrino mass spectrum: 

where Xij = Mf/Mf, H = m} D m D and [16] 

f(x) = y/x (l - (l + a) ln(l + -) + —?—). (42) 
[ x 1 — X J 

CP asymmetry parameters £j are related to the leptonic asymmetry parameters through Yi 
as [UrSS] 

y L ^-^ = ^fgj (43) 

where is the lepton number density, ul is the anti-lepton number density, s is the entropy 
density, Ki is the dilution factor for the CP asymmetry e$ and is the effective number 
of degrees of freedom |50j at temperature T = Mj. The baryon asymmetry Yg produced 
through the sphaleron transmutation of Yl, while the quantum number B — L remains 
conserved, is given by [HI] 

Y B = - *Jl±^E Yl (44 ) 
where Np is the number of fermion families and Nh is the number of Higgs doublets. The 



quantity Y B = — §§Yl in eq. (44) for SM. Now we introduce the relation between Y B and i] B , 



where r\ B is the baryon number density over photon number density n 7 . The relation is 



Vb 



s 1 



1U 



Y B = 7.0394Yb, (45) 



o 
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where the zero indicates present time. Finally we have relation between tjb and £j 



T] B = -2.495 x 



9*i 



(46) 



This dilution factor K{ approximately given by [53H55 



1 8.25 (Ki 

— — h 

Ki K V 0.2 



1.16 



with Ki 



Hi 



(47) 



where Tj is the decay width of iVj and iJ» is Hubble constant at T = Mj. Their expressions 
are 

.M? 



huMi 



and Hi = 



1 4irv 2 """"^Mp 
where v = 246GeV and M P = 1.22 x 10 19 GeV. Thus we have 

M P H U 



(48) 



K 



1.66 x 4:^J~g~v 2 Mi 



(49) 



B. Calculation of lepton and baryon asymmetry with broken scaling ansatz 



The matrix H = m D mD shown in eq. (41) is important to study leptogenesis. For six 



possible rriD with broken scaling ansatz by parameter e are given in eq. [2TJ They will 
generate following six possible H in three pairs: 



m M 2 



^ o \ 

l+p 2 (l + k 2 ) + 2p 2 k 2 e Ipqe^-^il + k 2 + k 2 e) 
^0 lpqe-^- e \l + k 2 + k 2 e) q 2 l 2 (l + k 2 ) J 



with I 



Ma 
M 2 



m M 3 



\ 





q 2 l 2 (l + k 2 ) lpqe-W- e \l + k 2 + k 2 e 

^0 lpqe l ^- e \l + k 2 + k 2 e) 1 + p 2 (l + k 2 ) + 2p 2 k 2 e J 



with I 



M2 
M 3 



(50) 
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moM 1 



( 1 +p 2 {l + k 2 ) + 2p 2 k 2 e lpqe^-°\\ + k 2 + k 2 e) ^ 



K lpqe-^- e \l + k 2 + k 2 e) g 2 / 2 (l + fc 2 ) , 



with I 



M 
Mi 



m M 3 



( q 2 l 2 (l + k 2 ) Ipqe-^-^il + k 2 + k 2 e)^ 



I IpqeW-^l + k 2 + k 2 e) 1 + p 2 (l + A; 2 ) + 2p 2 k 2 e , 

( l+p 2 {l + k 2 ) + 2p 2 k 2 e lpqeW- e \l + k 2 + k 2 e) o\ 
lpqe-W- e \l + k 2 + k 2 e) q 2 l 2 (l + k 2 ) 





V 



with I 



with I 



W 3 



(51) 



M 



m Q M 2 



( 



\ 



q 2 l 2 (l + k 2 ) lpqe-W- e \l + k 2 + k 2 e) o\ 

lpqeW- e \l + k 2 + k 2 e) 1 +p 2 (l + k 2 ) + 2p 2 k 2 e 
Oy 



with I 



Mi 
M~ 2 



(52) 



Parameters in above six possible H are already defined in eq.(|9]) and only I is defined here 
along with every H. Interesting features of the three pairs of H are that for every pair one 
generation of right handed neutrino decouples and also its decay width vanishes and hence 
could not take part in generation of lepton asymmetry For the first pair Ni decouples, 
for the 2nd pair N2 decouples and for the 3rd pair N3 decouples. Apart from this one 
more interesting point is that first matrix of every pair have similar expression in their non- 
zero diagonal and off-diagonal elements whereas the 2nd matrix of every pair have similar 
expressions. So, we don't need to study all the three pairs. We will only study the first pair. 

First generation of right handed neutrino N\ decay width is zero. Lepton asymmetry is 
generated through decay of N2 and N3 only contribute. Decay asymmetries 62 arid £3 for 



the first form of the first pair in eq. (]50|) 

1 



£2 



£3 



Im( ^ 3) /(M 3 2 /M 2 ) 



Anv 2 H 2 2 
1 lm(H? 



Anv 2 H~ 



^f(M 2 /M 2 3 ) 



M 2 m 

AllV 2 

M 2 m 



4:71V 2 



Ff(l 4 ) 
F'f(l/l 4 ) 



(53) 
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where I = J jjf- and 



rl 2 p 2 (l + A; 2 ) sin a 



l + /c 2 + 



l+p 2 (l + k 2 ) 



l+p 2 (l + k 2 ) 

F' = rp 2 (l + k 2 )sin a = + ^ (54) 

(l + fc 2 )yp 4 + r 2 - 2rp 2 cos a 
The definition of different parameters for different mjj are given in eq. [9] and also we have 
used g 2 e 2i ( /3_e ) = re ta — p 2 . The washout factors for 2nd and 3rd generation are 

K2 = Ififi M A PH ^2M = 913 - 7 (^) t 1 + ^ + k2 + ^ 

1.66 x 47r v /^v 2 M 2 VeV/ L 

1.66 x 47r J~g^,v 2 M ; 



^ = = 913.7 (^) (1 + ^^ + ^-2^ cos a. 

1.66 x ATt+fqliv* M 3 VeV/ 



(55) 

where we have used v = 246GeV, M P = 1.22 x 10 19 GeV and g H = 110.25 for SM with two 
right handed neutrinos. With this washout factors we can determine the dilution factors «2 



and K3 using the formula given in eq. [47J Well equipped with the above formulae for 62, £3, 
K2 and K3 we can easily generate the expression for baryon asymmetry 

T) B = -2.495 x J2 — 

i 9 " 

= -2.27 x 1(T 2 \s 2 k 2 + e 3 K 3 ] . (56) 
An additional beauty is that the expressions for i]b for two matrices in a pair are same. 



For the 2nd matrix of the first pair in eq. (50), expressions for e 2 and K 2 are same as the 



expressions of £3 and K 3 for the first matrix of the pair and expressions for e 3 and K 3 are 
same as the expressions of e 2 and K 2 of first matrix of the pair. So, effectively t\b expression 
remains same. Consequence is same as for the first matrix of the pair. 

The expression of t]b depends on m , k, p, r, a, e and additional two parameters M 2 
and I = yj On the top of constrained parameter space from neutrino data, we have also 
explored the parameter space with the additional constraint arises due to baryon asymmetry 
5.5 x 10~ 10 < r] b < 7 x 10" 10 PHSH] for 0.1 < I < 0.9 and 1.1 < I < 10 (avoiding point 
of degeneracy / = 1) and 10 12 GeV < M 2 < 10 15 GeV. We have seen that change in the 
parameter space is negligible. Only sign of a is constrained for different I. For 1.1 < I < 1.54 
and 0.1 < / < 0.9 sign of a is negative and for 1.54 < I < 10 sign of a is positive. Again 
a = 0, ±180 are not allowed. But value of a near and 180 are still allowed and the M 2 
value is large there, M 2 ~ O(10 15 ) GeV. 
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VII. SUMMARY 



To sum up, we have explored a predictive and testable scenario of neutrino mass matrix 
accommodating scaling ansatz with four zero Yukawa textures advocating type I seesaw 
mechanism with diagonal charged leptons and right chiral neutrino mass matrices. We 
break scaling ansatz in the Yukawa matrices to generate nonzero #13 through a dimensionless 
parameter e. The parameter space of the textures studied allow the 3<r value of 6*13 along 
with other neutrino experimental data. Using the #13 constraint we have restricted Dirac 
CP phase 5 and Jqp- We have also studied baryogenesis via leptogenesis arises in those 
textures, however, there is no drastic change in the parameter space due to the constraint 
from baryogenesis. But, the sign of the only phase present in this model is fixed. 
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Appendix A: RG Effect 



It is to be noted that even after breaking of scaling anasatz, the textures given in eqn.(21 ) 
are invariant under RG evolution. This is guaranteed in the following way : 

Following the methodology presented in Ref. [59]- [60] due to r lepton mass correction on 
m D oi eqn. ^ with scaling ansatz we get 



1 
1 
1 — A T 



(X\ CL 2 ^3 

kbi kb 2 kb 3 
h b 2 b. 



\ 



( 



\ 



d\ CL 2 

kb\ kb 2 kb$ 

) \0i b 2 o 3 J A r) 62(1 -A T ) 63(1- A T )j 

Redifining 6/s as &i(l - A T ) ->■ b x , b 2 (l - A T ) ->■ b 2 , 6 3 (1 - A T ) ->■ 63 we get 

( n n \ 



(Al) 



a x a 2 a 3 

kh(l + A T ) kb 2 (l + A T ) kb 3 (l + A 7 
h b 2 63 



(A2) 



/ 



where we consider (1 — A T ) 1 ~ 1 + A T since A T is far less than unity. If we consider 
k(l + A r ) — > k then, we get the structure of mo given in eqn.(|5]). So, m D with scaling 
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ansatz remains form invariant including RG effect. 



Now, if we consider scaling ansatz breaking through e parameter, the structure of comes 
out as 



' cii a 2 a 3 » 



kbi kb 2 (l + e) kb 3 
h b 2 b 3 

Again, RG effect through parameter A T on tud with broken scaling ansatz is given by 



V 



(A3) 



/ 



10 » 
1 
(1- A) J 



( 



d\ (X2 ^3 

kbi kb 2 {l + e) kb 3 
h b 2 b 3 



\ 



at a 2 
kh kb 2 (l + e) 
yhil-Ar) b 2 {l-A T ) b 3 {l-A T )J 



«3 

kb 3 



\ 



■ (A4) 



Performing the same exercise of redefinition of frj's and k, the same is obtained as in 



eq.(A3). So, the structure of mo matrices with broken scaling ansatz in eq.(21) are free 



from RG effect. 



Appendix B: Expressions used in Sec-4 



In our calculation we have written m v by breaking it into two parts, i.e 



m u = m° + em' u . 



(Bl) 



If we assume a generic form of m' v as 



( A x B 1 C x ^ 



B\ B 2 C 2 
Cx C 2 C3 



(B2) 
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(In our case A x = 0, B x = kp, d = 0, B 2 = 2k 2 p 2 , C 2 = kp 2 , C 3 = 0.) 

The different elements of the matrix h p in terms of the parameters (k, p, r, a) are given by 

h p n = rr^{2Re{A 1 ) + 2kpRe(B 1 ) + 2pRe(Cl)) (B3) 

h\ 2 = m 2 (Bl + kpB* + pC 2 + A x kp + B x k 2 re~ ia + dkre'^) (B4) 

h p 3 = ml{C{ + kpC* 2 + pC* 3 + A lP + B x kre~ ia + C x re~ i(x ) (B5) 
h p 22 = ml{2kpRe{B 1 ) + 2k 2 r(Re(B 2 ) cos a + im(-B 2 ) sin a) + 2kr(Re(C 2 ) cos a 

+im(C 2 ) sin a)) (B6) 

/4 = m 2 {kpC{ + A; 2 re* a C* + A;re ia C* + £ lP + B 2 kre~ ia + C 2 re- ia ) (B7) 
ft| 3 = ml(2pRe(Ci) + 2kr(Re(C 2 ) cos a + Im(C 2 ) sin a) + 2r( J Re(C 3 ) cos a 

+im(C 3 ) sin a)) (B8) 

Parameters like x, y, z, etc can be expressed in terms of different elements of h p matrix as 

x = c 12 (h^c 12 - ^ 2 e^c 23 si2 + ^ 3 e l0 s 23 si 2 ) - c 23 si 2 (/i? 2 e~^ci 2 - c 23 si 2 /i| 2 + hl 3 s 23 s 12 ) 
+s 23 s 12 (h( 3 e- i ' t 'c 12 - h( 3 e-^c 23 s 12 + h p 33 s 23 s 12 ) (B9) 
y = c 12 (h p u s 12 + h p l2 e l(t> ci 2 c 23 - h p l3 e l<t> s 23 c 12 ) - c 23 s 12 (h p 12 e~ l(l) s 12 + h 22 c 12 c 23 - h p 23 s 23 c 12 ) 
+s 23 s l2 (h{ 3 e-^s 12 + h p le~^c 12 c 23 - h p 33 s 23 c 12 ) (BIO) 
z = c 12 (h p 12 e l4 's 23 + h p 13 e t<p c 23 ) - c 23 s 12 (h^ 2 s 23 + h p 23 c 23 ) + s 23 s 12 (h p 23 e~ l4 's 23 + h p 33 c 23 )(Bll) 
w = s 12 (h p n si 2 + h p 12 e l(l> c 12 c 23 - h p 13 e t<p s 23 c 12 ) + c 12 c 23 (h[ 2 e~ t4> s 12 + h 22 c 12 c 23 - h p 23 s 23 c 12 ) 
s 23 c 12 {h{ 3 e-^s 12 + h p le-^c 12 c 23 - h p 33 s 23 c l2 ) (B12) 
q = si 2 (h p l2 e %<t> s 23 + h p l3 e l(l> c 23 ) + ci 2 c 23 (/i 22 s 23 + /i 23 c 23 ) - s 2 3Ci 2 (/i 23 e~ M/> s 23 + ft| 3 c 23 )(B13) 

(B14) 
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